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The motion of a small compact (SCO) object in a background spacetime is investigated further 

in the context of a class of model nonlinear scalar field theories that have a perturbative structure 

analogous to the General Relativistic description of extreme mass ratio inspirals (EMRIs). We derive 

regular expressions for the scalar perturbations generated by the motion of the compact object that 

are valid through third order in e, the size of the SCO to the background curvature length scale. 

Our results for the field perturbations are compared to those calculated through second order in e by 

Rosenthal in ill and found to agree. However, our procedure for regularizing the scalar perturbations 

is considerably simpler. Following the Detweiler-Whiting scheme, we use our results for the regular 

expressions for the field and derive the regular self- force corrections through third order. We find 

agreement with our previous derivation based on a variational principle of an effective action for the 

worldline associated with the SCO thereby demonstrating the internal consistency of our formalism. 

This also explicitly demonstrates that the Detweiler-Whiting decomposition of Green's functions 

is a valid and practical method of self force computation at higher orders in perturbation theory 

'■^ and, more generally, at all orders in perturbation theory, as we show in an appendix. Finally, we 

I ^ identify a central quantity, which we call a master source, from which all other physically relevant 

quantities are derivable. Specifically, knowing the master source through some order in e allows one 

^i to construct the waveform measured by an observer, the regular part of the field and its derivative 

^_^ on the worldline, the regular part of the self force, and various orbital quantities such as shifts of 

fi the innermost stable circular orbit, etc., when restricting to conservative dynamics. The existence 

^^ of a master source together with the regularization methods implemented in this series should be 

I indispensable for derivations of higher-order gravitational self force corrections in the future. 

I. INTRODUCTION 

> . . 

\^ This paper continues a study begun in |2] on higher-order self force effects in a nonlinear scalar model as an analog 

\^ of extreme mass ratio inspirals (EMRIs). In [2], henceforth referred to as Paper 1, we derive the finite expressions 

r — ■ for the self force on a small compact object (SCO) moving in a background spacetime (such as a supermassive black 

^^ hole) through third order in the small ratio e = Rm/Ti-- Rm is the typical scale associated with the size of the SCO 

f — . and TZ is the curvature length scale of the background spacetime. If the SCO with mass m moves in the strong field 

f^ region of a supermassive black hole with mass M then e ~ m/M is the mass ratio of the two bodies. 

^~~l There are several reasons for studying higher-order self force effects |2] ■ These include a need for higher-order self 

'""' force corrections to maintain the consistency of the inspiral's adiabatic evolution [2j and a need for high-accuracy 

^ waveforms to accurately measure the parameters of detected EMRIs in order to perform precise tests of General 

Relativity with space-based gravitational wave detectors. In addition, the recent discovery of transient resonances, 

which can significantly dephase the gravitational waveforms if not taken into account, requires modeling EMRI sources 

jrt with self force corrections through at least second order in e 0]. Higher order self force corrections might also be 

useful to describe binaries with less than extreme mass ratios, such as binaries with intermediate mass ratios in the 

range of 10" ^ - 10""*. 

In this paper we continue to focus on higher-order self force effects in a nonlinear scalar model of EMRIs [5]. 
Historically, (linear) scalar models offer a simpler framework for studying underlying issues of self force regularization 
and for developing practical self force computational schemes. Because of the relative simplicities that scalar models 
afford, it seems likely that the physics of higher-order self force effects can be investigated more easily and more 
quickly than in the gravitational EMRI context. 

In addition to the reasons given above and in Paper 1 for studying higher-order self force effects, our motivation 
for this second paper comes from several other directions. The first touches on issues of regularity for self force com- 
putations. In 5j, Detweiler and Whiting introduce a physically well- motivated decomposition of the retarded Green's 
function into regular and singular parts. The regular part satisfies a homogeneous wave equation and is responsible for 
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generating the regular part of the self force on the SCO. The singular part satisfies an inhomogeneous wave equation, 
carries all the information about the singular structure (that arises from using a point particle representation for the 
SCO) and exerts no force on the SCO. 

This decomposition into regular and singular parts is particularly useful for several reasons, including the property 
that the regular part of the field is differentiable on the worldline. Consequently, one needs only to compute the 
regular part of the field and its derivative on the worldline in order to evaluate the regular part of the self force. This 
approach to the self force problem, which we refer to as the "Detweiler- Whiting (DW) scheme" , is demonstrated in ^5] 
only at first order in e. The question naturally arises as to whether or not higher-order self force corrections can be 
computed in a similar manner. We address this question here and explicitly show that the DW scheme for calculating 
the regular part of the self force holds through third order in e and is, in fact, a well-defined and valid procedure at 
all orders in perturbation theory, at least in the class of nonlinear scalar models considered here. 

The second motivation comes from a desire to simplify the regularization and renormalization procedures that tend 
to impede progress in (higher-order) perturbative calculations of the self force. To our knowledge, Rosenthal is the 
only one to derive regular expressions for the nonlinear scalar field perturbations generated by the motion of a point 
particle in a curved, background spacetime through second order in e [Tj. As first pointed out using a specific nonlinear 
field theory in [T], scalar perturbations at second order in e involve integrals that diverge everywhere in spacetime. 
Rosenthal's method to obtain well-behaved, finite solutions to the field equation is based on a rather complicated and 
technical series of steps, which are not so easily generalized to gravitational perturbations because of gauge issues [6] 
or to apply at higher-orders in e. 

The regularization and renormalization methods we use here and in Paper 1 are borrowed from quantum field 
theory j7j , which has a long history of handling divergences and making sensible, finite predictions that are routinely 
confirmed by experimental measurements. These techniques are also applicable in curved spacetimes [5]. The theory 
is well-developed so that there are several practical and effective schemes for regularizing divergent integrals, one 
of which is dimensional regularization |^. Dimensional regularization takes advantage of the observation that the 
degree of singularity of divergent integrals depends on the number of space-time dimensions, d. In some dimensions 
the integral is actually convergent, which allows one to analytically continue the resulting finite expression for the 
evaluated integral to the physical number of dimensions^. 

We will not use dimensional regularization in most of our calculations in this paper but instead will explicitly 
carry around the singular integrals to demonstrate that these divergences can be cancelled in a consistent manner 
by introducing counter terms into the action [71 [5]. If, instead, we regulate these singular integrals in dimensional 
regularization then, as shown in Paper 1, these divergences have a vanishing finite part as a result of the analytic 
continuation implicit in this method. However, even without utilizing dimensional regularization, our approach is 
considerably simpler and more practical than that of [1]. 

The third motivation is a technical one. The methods of Paper 1 may be unfamiliar to researchers experienced with 
traditional methods of perturbation theory. Hence, in this paper, we approach the self force problem following the 
DW scheme. We also wish to show that our formalism is internally consistent and derive the regular part of the self 
force through third-order in e directly from the regular part of the field instead of from a variational principle for the 
effective action, as in Paper 1. A consistency check of the formalism is that the counter terms derived here are the 
same as in Paper 1. Indeed, we find precisely the same counter terms and exactly the same expression for the regular 
self force through third order. The DW scheme is a more direct and preferable way to derive self force expressions, 
in our opinion, than from the effective action of Paper 1. Consequently, we hope the methods used in this paper are 
more familiar and useful than those in [2]. 

Our last motivation for this second paper comes from what we think is an interesting theoretical observation. We 
are able to identify what we call a master source S that appears as the central, fundamental quantity from which 
most (if not all) other physically relevant quantities are derivable. In particular, knowing the master source through 
some order in e allows one to calculate the (scalar) waveform that an observer would measure in their detector, the 
regular part of the field anywhere in spacetime (including on the worldline), and the regular part of the self force. In 
addition, if one restricts oneself to a conservative description for the system then the corresponding master source can 
be used to compute shifts of the inner-most stable circular orbit, the energy, angular momentum, orbital frequency, 
and the redshift factor u* associated with the SCO's orbit, etc. The master source is straightforward to calculate 
using the methods of this paper. As a demonstration, we derive the expression for the master source through fourth 
order in e, which one then can use to find regular expressions for the field, the waveform, the self force, etc. 



^ It is well-known that there may appear singularities when continuing back to 3 -(- 1 dimensions. However, these divergences always 
manifest as poles in the complex plane of d, which can be absorbed into counter terms. In this paper, the analytic continuation of 
integrals that nominally diverge as a positive power of some cut-off regulator gives zero in 3 -f 1 dimensions. See [2] [8] for details. 



A. Organization 



This paper is organized as follows. In Section I B we review the nonlinear scalar model of EMRIs, which is meant to 
serve as an analog of the corresponding gravitational system. In Section [l C| we highlight the main results from Paper 
1 [2]- In Section Owe show how scalar perturbations are computed in the effective field theory (EFT) approach using 



Feynman diagrams. In Section III we compute the scalar perturbations emitted by the SCO through third order in 
e. We also regularize the singular integrals and renormalize the theory by introducing counter terms into the action. 
In Section IV we compare our results through second order in e to those of Rosenthal ^ . In Section W\ we derive the 



regular part of the self force in the DW scheme, i.e., from the regular part of the scalar field perturbations previously 



calculated in Section III We find exactly the same self force expression that we found in Paper 1 thereby providing 
direct evidence for the consistency of our formalism. In Section |VI| we identify the master source and calculate it 
through fourth order in e. 

We work in units where G = c = 1, in which case the small parameter for perturbation theory is e ~ m/TZ if 
m ~ Rm is the mass of the small compact object. Throughout, we ignore finite-size effects (such as tidal moments 
induced on the SCO) and spin angular momentum for the SCO, which will be considered in future work. Notation is 
as explained in Paper 1. 

B. The nonlinear scalar model of EMRIs 

In Paper I we introduced a class of nonlinear scalar theories designed to have the same structure as the perturbative 
description of EMRIs in General Relativity. Specifically, if the gravitational waves are in the Lorenz gauge and the 
SCO moves in a vacuum background spacetime, then the nonlinear scalar model that most resembles the perturbative 
structure of gravitational EMRIs has an action given by [5] 

S[z^^A]^-\J4>,o.<i>'''A\ci,)-mjdTB{4>) (1.1) 

Here, J = / d'^xg'^/'^{x)^ J dr = J^ dr, and A^ and B have a series representation in (jj 
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^('^) = 1 + E;^'^" (1-3) 

n=l 



where {a„,&„} are freely specifiable coupling constants or parameters. The equations of motion derived from (1.1) 
are 

a^ = -P^''V^lnB(0) (1.5) 

where a prime on A or _B denotes differentiation with respect to (j). 
Interestingly, defining a new field variable tp that satisfies 

V'.a = 0,a^(</') (1.6) 

leads to a field theory that is linear in the sense that there are no field self-interactions since such terms proportional 
to V'" for n > 2 are absent from the ensuing action, 

S[z'',ij]^-^J^,^i^'"-mJdTC{tP) (1.7) 

Here, 

OO 

C 



c{^p) = i + Y,^r (1.8) 



n—1 



and the c„ can be expressed in terms of the {a„, 6„}. The equations of motion for the new field ip and the worldline 



follow from (1.7), 



U^ = mj dT^-^^^^^^p^C'{^) (1.9) 

a^ = -P^''V^lnC(V') (1.10) 

where a prime on C denotes differentiation with respect to ■0- 

C. Selected results from Paper I 



In Paper 1, we introduced a new action principle, based on (1.7 1, that consistently incorporated the outgoing 
boundary conditions imposed on solutions to the field equation when integrating out the field to obtain the effective 
action. The variation of the resulting effective action yielded the worldline equations of motion and hence the self 
force to the given order in e. Specifically, expressions for the self force through third order in e were derived in [2] 
(in terms of the retarded propagator, which is generally unknown in closed form for most spacetimes) and found to 
contain singular contributions that depended on the worldline's past history. These divergences were shown to vanish 
in dimensional regularization. However, the formal expressions for those divergences were carried explicitly through 
the calculations to demonstrate that they could be absorbed into counter terms for the appropriate parameters of the 



theory. Specifically, these singular pieces were removed by adding the following three local counter terms to (1.7), 



dT-Si f dTiPiz'')-^ f dT^/{z'') (1.11) 

where 

,_^..(A),.3,^,(A)%^(A)^ (1.13) 

^..^^c^(A)+,..,,e3(A) (1.14) 

Then, the remaining regular part of the self force (using the Detweiler- Whiting decomposition for the retarded Green's 
function [5]) was found to be 

F]^(t) = (a^ + P^'^V,)|m2c2/;,.(z^') - ™'c2c2 Q4(z^) + y dT'i^K(^^2^')/i^(^^')') 

+ m^cl4 (inizn J dr' Dn{z^, z^^')Ir{z^') + j dr' dr" Dr{z^^, z^^')Dr{z^^' , z^")In{z'^")] 

Setting F^ equal to ma^ yielded the equations of motion for the SCO through third order in e. 

We also showed in Paper 1 that collecting all the terms proportional to the acceleration together gave a compact 
representation for the worldline equations of motion (i.e., self force) 

mr(z'')a^ = -mP^''V^r(z^) (1.16) 

where 

r(z^) = 1 - mcllniz'^) + m^clc2 Q^X^^) + J dr' i^^(z^ z'^')/^(z^V 



m^c^-^ 



IcI(ir{z^) j dr' Dr{z^,z^^')Ir{z^') + j dr'dr'' DR{z^',z^')DR{z^^\z^'')Iri{z^''')\ 



^^ (^4(-") + / dr' DRiz^, z^')ll{z^') + 0{e*) (1.17) 



Equivalently, dividing both sides of (1.16) by mr{z'^) gave 



(1.18) 



so that a comparison with ( 1.10 ) imphed T{z'^) = C{iPr{z'^)) that, if true, suggested that the regular part of the field 
would be given by 



(1.19) 



i^R{x)^ -mcilR{x) + m'ciC2 / dr Dnix, z'')Ir{z^') - m-'cici / drdr' Dnix, z^')DR{zf', z^ )lR{zf' ) 
""'"^'^^ I dTDR{x,zniU^n + 0{e') 



We will show in this paper that T{z'^) = C {tf) b.{z^^)) and that (1.19) is indeed the correct expression for the regular 
part of the field through third order in e. 



II. SCALAR PERTURBATIONS IN THE EFT APPROACH 



In this section, we describe how scalar field perturbations are computed within the EFT framework via diagrammatic 
methods. Perturbations of the scalar field ip{x) are the perturbative solutions to the wave equation ( 1.9 ) with a source 
due to the coupling between the field and the worldline for the SCO 



□i/) — m dr 



S'^ix^' - z^'{T)) 



jl/2 



1 2 

ci + C21JJ + -csijj + 



(2.1) 



If D{x,x') is the Green's function, or "propagator", appropriate to the problem (e.g., for outgoing boundary condi- 
tions) then the solution to (2.1) is formally 



^(a;) = -mci dr D{x,z'') - mc2 dr D{x, z'')^P{z^') - -mc^ dr D{x, z^')^p^iz^') + ■ ■ ■ (2.2) 

The perturbative (but divergent) solution can be found by solving the above expression iteratively, which through 
second order in e is 



i:{x)^~mci / dTD{x,z'') + m^ciC2 / drdr'Dix, z'')D{z'' , zf"') + 0{e^) 



(2.3) 



with the higher-order terms becoming increasingly more complicated and more involved to compute^. In addition, 
one has to carry around many terms just to find the 0(e") contribution for some n. For these reasons, it is beneficial 
to perturbatively solve the wave equation (2.1 ) using Feynman diagrams, which allow one to systematically compute 



only those contributions appearing at a given order without the need to carry around many terms that may otherwise 
be irrelevant at that order. 

Feynman diagrams are constructed by drawing all connected, tree (i.e., not containing any loops of propagator lines, 
which correspond to quantum corrections) diagrams that scale with e to the given power. A diagram representing 
a field perturbation must have one "dangling" propagator (i.e., curly) line having one end that is not attached to a 
worldline. The free end of that line represents the field point where one is measuring the field. To determine the 
scaling of a diagram, which is a procedure called power counting, we recall that the scalar perturbations vary with 
the scale of the background spacetime curvature so that datp/^p ^ TZ~^ and hence a;'" ^ TZ. In units where G = c= 1, 
the action has dimensions of (Length)^ from which it follows that ip is dimensionless in 3 -1- 1 spacetime dimensions 
and ^ ~ 7?,°. The retarded propagator Drct{x,x') satisfies 



J X 



(2.4) 



^ It is not really the case here that the scalar perturbations are difficult to calculate through third order in e. However, calculating 
perturbative solutions for metric perturbations to some order in e quickly becomes nontrivial and quite involved, even for relatively low 
orders. We remind the reader that one purpose of this nonlinear scalar model is to provide useful insights into calculating higher-order 
self force corrections in the full gravitational EMRI case. 



FIG. 1: The Feynman diagram contributing to the radiated scalar perturbations at first order in e. This diagram scales as 
{Ti-~^){m){TZ) = £ where the m comes from the vertex, the propagator line scales as 1/7?.^, and the proper time integration 
scales as J dr ~ 7?.. 



from which it follows that D^ct{x,x') ~ 7?."^. Lastly, every worldline vertex in a diagram scales as '^ m (see first of 
the Feynman rules listed in AppendixfAl and / dr ~ 7?, is the typical dynamical time-scale for the motion of the SCO 
in the background spacetime. These scaling laws allow one to determine which order in e that a Feynman diagram 
appears in perturbation theory. An example applying the Feynman rules and power counting to a specific Feynman 
diagram is given in Appendix \K\ 

The fact that the perturbative solution for i^{x) in ( 2.3 1 can be expressed via Feynman diagrams is not so surprising. 
(This is well-known from quantum field theory.) In fact, (2.3 1 can be written more suggestively as 



il^{x)= / dTD{x,z^'){-mcl)+ i dr D{x, zf'){-mc2) / dr' L>(2^, z^ )(-mci) + ©(e^) (2.5) 

which is exactly the expression found by applying'^ the Feynman rules to Figures [I] and [2] ~ see (3.2) and (3.4| below. 



III. SCALAR PERTURBATIONS IN THE NONLINEAR MODEL 

In this section we calculate the regular part of the scalar perturbations emitted by the SCO through third order in e 
within the nonlinear scalar model. We show that all divergences can be removed by introducing suitable counter terms 
in the action. These counter terms will be shown to precisely equal those of Paper 1 despite the different contexts 
in which they appear. This thereby demonstrates the internal consistency of our regularization and renormalization 
program. 



A. Scalar perturbations at first order 

The leading order contribution to the radiated scalar perturbations comes from the diagram shown in Figure IT] 
According to the power counting rules this diagram scales linearly with e and yields the first-order contribution to 
the field. The Feynman rules imply that Figure IT] is equal to 

(3.1) 

(3.2) 

where D ^{x,x') = Diet{x,x'). Since we are assuming that the field point x^ is not on the worldline then the above 
integral is manifestly finite and describes the leading order scalar perturbations generated by the motion of the SCO 
as it moves through the background spacetime. 



V'(i)(a;) = I drD +(.t, z^)(-rnci) 
-mci / dr' D,:ct{x,z^ ) 



^ Start from the top of the propagator line in each figure, which is at the field point x'^, and worfc through the diagram, daisy-chaining 
propagators via worldline vertices and introducing the corresponding factors as stated in the Feynman rules in Appendix [A| 
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FIG. 2: The Feynman diagram contributing to the radiated scalar perturbations at second order in e. This diagram scales as 
(7^-2)(m)(7^-2)(m)(7^)2 = £^ 



B. Scalar perturbations at second order 

The diagram in Figure [2] scales as e^ and describes the scalar perturbations at second order. The Feynman rules 
imply that 

i^(2){x)= I dT'dr" D-+{x,z^''){-mc2)D-+{7y ,zt'"){-mci) (3.3) 

= m^cicz / dT'dr" D,,t{x, z^'')D,,t{z^'' ,zf'") (3.4) 



(Compare (3.4) with the second order term in (2.3 1, which is found in a more traditional way.) The integral over t" 
diverges when r" = r'. 

To deal with this singularity, we use the Detweiler- Whiting prescription |S] where the retarded Green's function is 
decomposed into a regular (R) and singular (S) part, 



Drct{x, x') = Dr{x, x') + Ds{x, x') 



(3.5) 



While there are many ways to isolate the singular part from the regular one, this particular decomposition has the 
advantage that the regular part Djf satisfies the homogeneous wave equation and is the piece of the retarded Green's 
function actually responsible for describing the force on the SCO. However, the singular part Dg carries all the 
divergent structure of I?rot, satisfies the inhomogeneous wave equation and exerts no force on the SCO. 

In Paper I, we showed how to use the Detweiler- Whiting decomposition to isolate the singular part(s) of the worldline 
integrals encountered there. We refer the reader to Paper I for the details. To evaluate the divergent integral in (3.4) 
we recall from Paper I that 



where 



is the singular integral and 
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is the regular part of ( |3.6[ ). Then, ( |3.4[ ) equals 

dr'Dret{x,z''') + m^CiC2 I dr' D,,t{x, z''' )Ir{z''' ) 



'0(2) = m CiC2 
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(3.6) 



(3.7) 
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FIG. 3; The diagrams contributing to the radiated scalar perturbations at third order in e. Both diagrams scale as 



C. Scalar perturbations at third order 



At third order in e there are two diagrams contributing to the emitted field, which are shown in Figure [3] The 
Feynman rules imply that 



^(3)(a;)= dT'dT"dT"' D-+{x,z>'){^mc2)D-+{zi' .z^" )(-mc2)L'"+(z'" , z^ )(-mci) 



+ - / dT'dT"dT"' D-+{x,z^' ){~mc3)D-+{z^' ,2^ ){-mci)D-+{z'' , z^ )(-mci) (3.10) 

where the factor of 1/2! in the second term comes from the symmetry factor of the diagram on the right in Figure p] 
With D-+{x,x') = D^etix,x') it follows that 

,3„ 2 



^(3)(.x)= -m'ci4 I dT'D,,t{x,z^) I dT"dT"'D,,t{z\z^' )A-ct(^^ ,z^ ) 

-[ 2 



1 



m C1C3 dr' Drct{x,z'^ ) 



dT"Act(^^z^ ) 



(3.11) 



From Paper I we recall that 

dT"dT"' Act {z''' , z^'" )D,,t {z^" , 2^'" ) 



47r / 



I +2[^]ir{z^')+ I dT"Dn{z^'z^')In{z^^^^) (3.12) 



4tt J' 



Therefore, using ( 3.12| in the first line of (3.11) and (3.61 in the second line gives the following expression for the 
contribution to the field at third order, 



V-o) = -m^cicl I dr' D,^t{x,z''')\ f — ) +2 
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yR{z^')+ j 



dT"DR{z^ ,z^ )Ir{z^ ) 



m C1C3 / dr' D ret {x 



■ =''){{i) +2(A),,,„.', + 4„.^, 



(3.13) 



D. Renormalization 



Combining our results for the first, second and third order contributions to the scalar perturbations emitted by the 
SCO from (3.2 1, (3.4), and (3.131 gives the solution to (2.1 1 through third order in e. 



ip{x) — I dr' D,-ct{x, z^ )<^m C1C2 ( ";~ I ~ m CiC 
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(3.14) 
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FIG. 4: The contributions of the counter terms to the radiated scalar perturbations through third order in e. 



As with the self force calculation in Paper 1, we find singular contributions to the radiated field that are of higher- 
order (e.g., ^ A^) and history-dependent (e.g., ^ AIji{z^ )). Despite the appearance of these seemingly pathological 
structures we can cancel their contributions by adding purely local counter terms into the action. In particular, we 
add to the action (|1.7[) the following counter terms 



-(5i / dTV(2^)-y / drV'(z'^) 



(3.15) 



since these are the ones at lowest orders in e that describe a coupling between the worldline and the scalar perturbation. 
Consequently, the counter terms modify the equations of motion and their (divergent) solutions so as to make them 
finite. The contributions to the radiated scalar perturbations from these counter terms are displayed in Figure HI 

Following steps similar to those used in deriving (3.14) we find that the contribution of these counter terms to the 
field is 

V'Ct(2:) = dr' D^ctix,z^' )< ~ Si+mciS2i—j +mc2Sii — 

+ f dr' D,,tix,z^')lR{z''')l^mciS2+mc2Si\ +0{e'') (3.16) 

Requiring that ip + ipcT be regular and independent of the divergent quantities A implies that i5i and 62 are given by 

A\ c. Q / A \ m^clcs f A 
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(3.17) 
(3.18) 



Therefore, the regular scalar perturbations emitted by the SCO through third order in e are 

'(/'rad(a:^) = 'ipix) + tpcT{x) 

= I dr D,^t{x,z''')\ -mci+m^ciC2lR{z''')-m^cicl [ dr" Dr{z>'' , z''")Ir{z''") 



(3.19) 



m^clc3 2 



/A(^^ ) 



+ 0{e') 



(3.20) 



This is one of our main results. The field i/'rad in (3.20), which we will call the radiative field for definiteness, is the 



regular field measured by an observer away from the SCO and, by analogy, corresponds to the gravitational waveform 
of an EMRI source. It is a strong check of self-consistency that the counter term coefficients 61 and S2 in (3.17) and 



(3.18) have the same values as first derived in Paper I (and given in (1.131 and (1.14)) because of the different ways 



that they are obtained. Specifically, the counter terms here are introduced to make the field regular while those in 
Paper 1 to make the self force regular. 



For later convenience, we will write (3.201 in a more compact form 



V'rad(2:) ^ / dT'D,,t{x,z'' )Sr{z^ ) 



(3.21) 



where the regular source term Sji{z^ ) is 



.3^2, 



Sr{z'^) = -mci + m^ciC2lRizn - m^c^cl / dr' Dr{z^ ,z''')Ir{z^^') - "^|i^4(z'^) + 0{e^) (3.22) 
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IV. COMPARISON WITH ROSENTHAL'S SECOND ORDER SCALAR PERTURBATIONS 



Given that (3.20) is the expression for the regular part of the radiative field and is valid through 0{e^), we next 
compare our result with the O(e^) regular field derived by Rosenthal in [L . There, a particular nonlinear scalar field 
theory is chosen such that it yields the following field equation 



□ c 



dr 



(54(.t^-z^(t)) 
oV5 



(4.1) 



where q denotes the scalar charge of the particle. The challenge in solving this equation, as noted by Rosenthal [T] , 
can be seen by decomposing the field as 



(x) = #i(x) + q^Mx) + 0{q^) 



(4.2) 



where the subscript labels the order of the perturbative solution. Substituting this into (4.1) and using the solution 
for the first order perturbation 

(4.3) 



(4.4) 



01 (a;) = / drA-ot(a;,z^) 
gives the following equation for 02 

n02 = /"drdr'V„i?,et(a;,z^)V"A-ct(a;,z^') 
The corresponding solution is formally given by 



M^) = - / Act(x,y) / drdr' V^, Act (y,^^)V"Aet(2/, ^'^ ) 



(4.5) 



which diverges as the inverse second power of the proper spatial distance orthogonal to the worldline as one integrates 
over 2/'^ near the particle. In fact, this divergence arises at every field point a;^ in the spacetime. 

In [T] Rosenthal devised a series of steps to obtain well-defined solutions to (4.1 1. Through second order, his 
prescription yields 



(t>ix) 



q I dTD,,tix,z^) + ^ /dTAct(a:,^^) - <z' / 



drAct(a;,z^)/fl(z^) + 0(e') 



(4.6) 



where Ir{x) is given in (3.8 1. Rosenthal derived this solution by carefully following several steps: (1) investigating 



the behavior of the wave equation when q — > 0, (2) making an ansatz for the particular solution to the wave equation, 
(3) using physical considerations to identify the divergent boundary conditions for the field (/>(a;) as the field point 
x^ approaches the worldline, and (4) solving the wave equation with the divergent boundary conditions from the 
previous step. This is a rather tedious process in practice and becomes more involved when calculating second order 
gravitational perturbations [BJ or going to higher-orders in e. We show below that one need not employ such a series 
of steps. Instead, we will recover Rosenthal's perturbative solution for the field (4.6 1 using effective field theory 
techniques. 

Rosenthal's nonlinear field theory is a member of our class of nonlinear scalar models. To make a meaningful 
co mpar ison, wc need first to find the values of th e pa ramet ers { ambn} appearing in (1.4). From the wave equation 
in (1.4 1 and the definitions of A^{(f)) and i?(0) in (1.2 1 and (1.3), respectively, it is easy to show that 



~A 

l2 



6 

= h- 



h- 



a^bi 



(4.7) 
(4.8) 



Inserting (4.7 1 and (4.8) into the wave equation for (j){x) in ( |1.4[ ) and comparing the resulting expression order-by-order 
in powers of (j) with Rosenthal's wave equation (4.1) implies that 



bi = 

b2 = 
as = 



(4.9) 

(4.10) 
(4.11) 
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and a„>4 = 0, &n>3 — 0. Substituting these values for {a„, 6„} into the renormalized expression for the radiative field 
given in (3.20) gives 



^rad(a;)= / dTD,,t{x,z^)\q-q'lR{z>') + 0{e-') 



Using the expression for ^ in terms of ■0 = V'rad by applying the inverse of the field redefinition of (1.6) yields 



1 



(x) = V^radla:^) + ^rV'radCa;) + 



so that substituting in for VAadl^;) from (4.12) gives 



drDrctix^z'') 



drA-ct(a;,z^)/fl(z^) + 0(e3) 



(4.12) 



(4.13) 



(4.14) 



Comparing with (4.6) we see that we have reproduced exactly Rosenthal's result [T. However, we have arrived at 
(4.14) using a simpler and more straightforward regularization procedure than that developed in [T]. 

Our approach to regularization can be streamlined even further by evaluating the singular integrals in dimensional 
regularization. Since all such integrals are not logarithmically divergent then the divergences will vanish in dimensional 
regularization. Regularizing in this way then becomes rather trivial. 



V. SELF FORCE FROM BACKREACTION OF SCALAR PERTURBATIONS 



In the previous section, we found the formally divergent solution to the wave equation for iIj{x), given by (3.14), 
through third order in e. The divergences were removed by introducing counter terms into the action, which yielded 
the regular scalar perturbations radiated to an observer. In this section, we show how to obtain the regular part of 
the self force through third order from knowledge of ^rad(3;) (i.e., we implement the DW scheme through third order). 
We also compare with the result we obtained earlier in Paper 1 for the self force. We show in Appendix [B] that the 
DW scheme can be applied at any order in perturbation theory. 



A. Regular part of the self force 



In renormalizing the field ip{x) in Section III we added the counter terms in (3.15) to the action (1.7), 

S[z'',i;] + ScT[z^,tp] ^-\j V'.aV''" - / drlm + (mci + 61)4, + ^(mc2 + (Ja)^' + ^mcgV'^ + 0{e^)\ (5.1) 

In doing so the field and worldline equations of motion acquire additional contributions arising from the presence of 
these counter terms. Specifically, the wave equation becomes 



□V'(^) = / dr ^-^^^""^^p^Umc, + <5i) + [mc^ + 5,)^ + "^^^ ^ ^^^4) 



(5.2) 



The resulting solution to this wave equation is precisel y '4'{x ) — fAradC^), the radiative field already given in (3.20), 
which can be verified by direct substitution and using (3.17) and (3.18). However, this should be obvious since the 



counter terms are introduced in the first place to cancel the divergences in (3.14) and yield the finite expression in 



(3.20), namely, ■(/'rad(a;) 



The self force also inherits extra contributions from the counter terms. In addition, the self force that results by 
varying the action in (1.7) with respect to the worldline coordinates is in terms of the radiative field, as this is the 



quantity that pcrturbatively solves the wave equation with the counter terms included, (5.2). The resulting self force 
expression is 

F^r) - - (a^ + P^"^V,)j(mci +5i)V,.ad(^^) + ^(mc2 + <52)V;L(^^) + ^^rad(^'^)| + 0(£^) (5.3) 



upon expanding out C{ip) using (1.8). Substituting in (3.20) for i/'rad(z'^) into F^^{t) in (5.3) then gives 



F^{T)^F^{r) 



%^c\ f A\ wt'c\c2 f ^Y rn^c\c\ / A \ ^ m^cXc'i / A ^ ^ 



47r 



47r 



2 V47I-7 ^ 6 



47r 



+ 0(£^) 



(5.4) 
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where F^{t) is the regular part of the self force through third order, 

+ m^cjcl flniz^ J dr' Dr{z'', z'^')Ir{z^') + j dr'dr" Dniz^', z^'')Dr{z^'' , z''")Ir{z^^")\ 

+ '^^{^^lU^n+ j dr' Dr{z^,z^')II{z^')^ +0{e^)^ (5.5) 



There are remaining divergent contributions in ( |5.4| that are not cancelled by the 5i and 82 counter terms. However, 
these divergences are all proportional to the particle's acceleration i mply ing that a mass counterterm 5m is sufficient 
to cancel these pieces. Therefore, adding —6m J dr to the action in (5.1) yields 



F^{T)-6ma^^F^iT)+a^ 



so that choosing 



(5m = 



- Sm. + 



^) + 



A 
in 



'm?c\c2 f A 
2 \4^ 



-l'--2 



A \ ^ „,4^3, / A \ 3 

47r 



m cfc^ / A 
6 V4^ 



0{e' 



n3^2^ 



?C2 / A \ 



An 



-V^'2 



A 
An 



m c\cz ( A 



6 



An 



+ 0(e4) 



renders the self force finite 



F^(r) - 5ma^ = F]^(r) 



(5.6) 
(5.7) 
(5.8) 



where F'j^ij) is given in (5.51 and is precisely the same expression derived in Paper 1 (shown in this paper in (1.15)). 
This is our second main result. In addition, the expression for the mass counter term in (5.7) is exactly the same as 



given in Paper 1. The agreement among the expressions for all three counter terms and for the regular part of the 
self force as derived here (via the scalar field perturbations) and in Paper 1 (via a variational principle) provides a 
strong check of the consistency, results, and tools used in the effective field theory approach for deriving self force 
corrections in a nonlinear field theory coupled to a worldline. 



B. An efficient derivation of tiie tliird order regular self force 



The regular part of the self force in (5.5) can be derived in a more straightforward way. All singular contributions 



vanish in dimensional regularization (i.e., A — 0), as already noted multiple times here and shown in Paper 1. 
Therefore, the counter terms (Si, (^2, and (S„i in p. 17 ), (3.181, and (5.7) all vanish, the regular part of the field 
evaluated on the worldline is just 



V'(^^) - ^rad(^^) - ^r{z^) = / dr' Dr{z'', z^ )Sr{z'^ ), 



and the formally divergent expression in (5.3) for the self force becomes manifestly finite. Therefore, 



mw 



F^{t) ^ F^ir) ^ - m{a'^ + P^'^V.) [c.i^R + |^|j + |v|, + 0{e') 



(5.9) 



(5.10) 



where i/'K is given in (5.9) and (3.22), which is subsequently substituted into (5.10) to recover the same expression 



for the regular self force through third order given in (5.5) 



VI. THE MASTER SOURCE 



The field ?A(a;) in (3.14) can be written as if there is a source 5(z'') generating the scalar perturbations where 



.( 



A\ 



S{zn^m'c,C2[--\-m'c,4\--\ 



'\An I 



A\' 



An I 



"^^'^s) -«<-'{ -=-<s)-"'=*3(^) 



n3^2. 



mci + m2ciC2/fl(z^) - m^cic^ / dr' Dr{z^ ,z^' )Ir{z^^' ) ~ '-'^^ll{z^) + 0{e^) 



(6.1) 
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so that 



V'(x)= dT'D,,t{x,z^')S{z^') 



(6.2) 



This representation for the field holds to all orders in e, which is easily seen by solving formally the exact field equation 
in((L9l), 



^{x) = -m / dr A-ct(x,z'^)C'(V'(z'')) 



implying that 



5(z'^) = -mC"(V'(z^)) 



(6.3) 



(6.4) 



For reasons that will become clear below we call iS(z^) the master source function^. It will be convenient to identify 
the regular and singular parts of the master source, 5(z^) = Sr{z^^) + Ss{z'^), with 



5fl(z^) - -mC'ii^Riz^) 
Ssizn^Sizn-Sniz^) 



(6.5) 
(6.6) 

(6.7) 



Through third order in e, the regular part of the master source Sn.iz'^) is given in (3.22) and the singular part by 



Ss{z'') = m^ciC2 ( ^ ) - "^ cic: 



A\ 



n3„ J2 



A \ m^c\c^ ( ^\ 



An J 



4ttJ 



I - lR{z'-)t - 2m'ct4 ( f ) - raV,cJ^ ) I + Ole') 



An 



V 47r 



In dimensional regularization A = and the singular part of the master source vanishes, Ss — 0. We will ignore Ss 
in the remainder. 

The key point to note is that the master source is the fundamental quantity to calculate for EMRIs: it is used to 
construct the radiative scalar perturbations observed by detectors ( 3.21| , 



the regular part of the field (5.9 1 



'(/'rad(a;) = dTD,ctix,z^')SB.iz^'), 



^Rix)^ / dTDR{x,z^')SRiz^), 



and, subsequently, the regular self force (1.10) 



n=l 



(6.9) 



(6.10) 



(6.11) 



where the self force depends on Sr{z'^) through (6.10) evaluated on the worldline. Also, the physical field is 4>ra.d{x) 



and the corresponding regular part on the worldline is (J)r{z^)^ both of which can be calculated by inverting the field 
redefinition in ( 1.6 ) in terms of V'rad(a^) and V'fi(^^), respectively, which both depend on the master source. Therefore, 
by computing Sr{z^) to the desired order in e one can compute essentially all relevant physical quantities describing 
the EMRI. This is our third main result. In the next subsection we show how to calculate the regular master source 
using diagrammatic techniques. 



* We call S{x) a "master source" since it is representative of all the field- worldline interactions that source the field measurable by some 
observer. Normally, S might be called an "effective source" but this term is already used to describe the finite, windowed source term for 
computing regular perturbations in some approaches to numerical self force computations (e.g., |10jV S might also be called a "dressed 
worldline vertex" as in |lll . 
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Cl 
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C3 



Cl 



(c) 



FIG. 5: The diagrams contributing to the regular part of the master source Sr(z^) at (a) first order (~ Tie), (b) second order 
{^ TZe^), and (c) third order (~ T^-e"^) in e. The solid gray line emanating from a vertex represents the place to attach a Green's 
function such as Drct or Dr in order to compute, for example, the radiated field V'rad(2;) or the regular part of the field tpR{x), 
respectively. 

A. The master source from Feynman diagrams through third order 



In this section wc demonstrate hovi' to eompute the master source through third order in s. The master source can 
be calculated through third order in e via the Feynman diagrams given in Figure [5] - these are all of the tree- level, 
connected diagrams that scale as e'^ or lower. The solid gray line attached to a vertex represents the place to connect 
an "external" Green's function such as D^etix, z^) if one wants to compute the radiated scalar field ij^radix) or Db,{z^) 
if one wants to compute the regular part of the field on the worldline 4'r{^^)j for example. Compare Figure p] with 
Figures [1][3] 

The master source at first order is given in Figure [5 (a) [ and is simply the Feynman rule for the SCO coupling to a 
single field perturbation. 



The master source has a contribution at second order given by Figure |5(b)[ which equals 

4'^(z^) = (-mc2) /"dr'Z?-+(z^,z^')(-mci), 

and at third order by Figure [5(c) 



(6.12) 

(6.13) 
(6.14) 



;^^\z^') = (-mc2) / dT'D-+{z^',z'''){-mc2)D-+{z^'\z^'"){-mcl) 

+ ^(-mc3) f dr' D'~+{z'',z^''){-mci) f dr" D-+{z^', z''"){~mci) 

3 2 

= — m C1C2 



j dT'DR{z^^,z'^')lR{z^^') - ^m'cjcslUzn 



(6.15) 
(6.16) 



In evaluating the diagrams in Figure p^ we use D '^{x,x') — Dfi{x,x') since all of the singular integrals vanish in 
dimensional regularization as we have pointed out multiple times. The master source through third order is then the 
sum of KT^, (pHI and KM, 



Sniz^ ^ -mci+ m^ciC2lR{zn ~ m^cicl / dr' i?fl(z^ ^^ )/fl(z'^ ) - -m^c^cg/Kz'^) + 0{e^), 



(6.17) 



which agrees with the expression we defined above in (3.22). From ( 6.17[ ) it is easy to obtain the radiated field tpra.d{x) 
in (3.20) using (6.9|, the regular part of the field on the worldline ^Pr(z^) in (5.9) using (6.IOI, and the regular part 



of the self force through third order in (5.5) by using (6.10 1 and (6.11). 



The master source at fourth order 



Perturbatively computing the master source is sufficiently straightforward that one may go on to compute the fourth 
order corrections and higher, if desired. Figure^ shows the four diagrams appearing at 0{£^). Using the Feynman 
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FIG. 6: The diagrams contributing to the regular part of the master source Sr{z'^) at fourth order in e. Each diagram scales 
as 7?.e*. 



rules to translate these diagrams into expressions gives 

Sr\z'') = {-mc2) [ dr' D-+{z'',z'''){-mc2) / rfr" Z3-+(z^', z^")(-mc2) j dr"' D-+{z^'" ,z^'"'){-mci) 



-A-mc2) I dr' D-+ {zt" , zt"' ){-mcz) ( dr" D-+ [z^"', z^"" ){-mcx) f dr'" D-+{z^'' , z''"'){~mci) 



2! 

+ (-mc3) / dr' D-+{z'',z^''){-mci) / dr" D-+{z^', z^"' ){-mc2) / dr'" D-+{z^'' , z^'" ){-mci) 

3 



+ ^{-mc4)(JdT'D-+{z'^,zy){-mc,)) 
or, again, with D ^(x,x') = Dii{x,x') in dimensional regularization, 

Sfiz^') = m'^ccl J dr'dr" Dr{z^,z^')Dr{z'^' ,z'^")In{z''") + ^m^cW^ j dr' Dr{z^ ,z^'' )ll{z^^') 



(6.18) 



(6.19) 



The field radiated to a far-away observer is then simply given by (6.9 ) and the regular part of the field on the worldline 



is (6.10) where the master source Sr through fourth order in e is the sum of the terms in (6.12), (6.14), (6.16) and 



(6.19). The contribution to the regular part of the self force at fourth order is found by reading off the appropriate 
terms from (6.11 ) upon substituting in (|6.10) in terms of the master source (we will not give the expression here). 



C. Discussion 

In the nonlinear scalar model considered here, the existence of a master source function is made possible, in part, 
because of the field redefinition from (j) to ip, which removed the self- interaction terms proportional to Df/i from the 
Lagrangian thereby ensuring that 



ip{x) 



dTD{x,znSR{z'') 



(6.20) 



for the appropriate propagator (e.g., retarded or i?-part). 

To see this, let us step back and consider the theory in the original variable 0(x). The action is given in (1.1) and 



the corresponding wave equation in (1.4). At second order in e there is a contribution from the diagram shown in 
Figure [7J which describe s no nlinear self- interaction of three scalar perturbations in the bulk spacetime. The following 
self-interaction term in (1.4) 



03 

'3! 



03 

'l2 



'ac 



(6.21) 
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FIG. 7: A contribution at second order to the perturbative solution of the wave equation for (f>{x). 



is the quantity that accounts for the interaction of the three propagator hnes in Figure [7] The following contribution 
to (/>(a;) at second order from Figure [t] is 

0(2)(a;) D / drdr' / Act(a;,2/) f y j D^Actly, 2'')(-m5i)i:'rct(y,2"')(-™6i) 



+ jj drdr' I UyD,,t{x,v)yj ]D,,tiy, z^){-mbi)D,,t{y, z''' ){-mh) 



(6.22) 



(We use D to indicate that the right side is only one part of the full expression for (^(2) (a;), the other part coming from 
evaluating a separate diagram.) Using the equations of motion for the retarded and advanced Green's functions, 



OyD,.ct{y,x') 



5\y-x') 

fll/2 



UyD^ct{x,y) = OyDa^dv{y,x) 



5\x ~ y) 

fll/2 



(6.23) 
(6.24) 



we see that (6.22) becomes 



(l^{2){x) D 



6 



drdr' D,,t{x,z^')DR{z^',z^') 
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drdr' D,,t{x,z^')D,,t(x,z^') 



(6.25) 



upon using dimensional regularization to evaluate the singular integral, which vanishes as usual. 

The key point is that the effect of D in (6.21) is to "collapse" a propagator line so that instead of the three 
propagators originally appearing in (6.22) one is left with only two propagators in (6.25). As a result, the last term 



of (6.25) contains a product of two Green's functions, each having one of its arguments evaluated at the field point 
x^ . Hence, the field 4){x) cannot be written at all orders in e simply as the integral of some source, as in (6.20). In 
fact, it is precisely the field redefinition from to ■0 in \1.&\ that eliminates those self-interaction terms proportional 
to Df/) (which are all of the self-interaction terms in this theory, in fact). Therefore, in the ij: variable there are no 



contributions like the last term in (6.25) and one can identify a master source 5, as in (6.20) 



Using another Green's function decomposition besides Detweiler and Whiting's, such as Hadamard's decomposition 
into direct and tail pieces, does not yield a master source (even after changing variable from to ip) since the tail 
field, while finite on the wordline, cannot simply be differentiated to yield the complete, regular self force (see, e.g., 
[2]). Consequently, the self force would not be expressable in terms of a master source. It seems to us that there is a 
special relationship between the field redefinition in (1.6) and the Detweiler- Whiting decomposition since they both 



conspire to reveal a master source Sr that generates many, if not all, physically relevant quantities for EMRIs. 



VII. CONCLUSION 



In this paper we derived the scalar perturbations generated by the motion of a small compact object in a class of 
nonlinear scalar models for extreme mass ratio inspirals that were developed and originally studied in Paper 1 of this 
series [5]. In particular, we showed here that one can obtain sensible and finite expressions for the field perturbations 
generated by the SCO through at least third order in e = m/TZ. We also showed that one can derive the self force 
corrections to the motion of the SCO from a knowledge of the regular part of the field perturbations when the latter 
are evaluated on the worldline. This shows, rather explicitly, that the Detweiler- Whiting scheme for computing self 
force corrections is valid through third order in perturbation theory. We also showed in Appendix [B] that the DW 



17 

scheme can be applied at any order in perturbation theory since the regular self force is found, at all orders in e, to 
be given by 

oo 

F^{r) = ~m(a^ + P'^'^V.) J^ SV'S(^^) (7-1) 

n. 

n—l 

The divergent integrals that we encountered in the course of the calculations here are subtracted by exactly the 
same counter terms found in Paper 1. This explicitly demonstrates the internal consistency of our regularization and 
renormalization program. If using dimensional regularization to regularize the integrals, then all singular contributions 
vanish since the nonlinear scalar model considered here cannot generate any logarithmically divergent integrals. This 



property greatly streamlines the derivation of higher-order corrections of the field and self force (see Section VB). 

For a specific choice of parameters in our nonlinear scalar model we recovered the second-order perturbative solution 
for the field (j){x) that was first derived using very different methods by Rosenthal [1]. The method of counter terms 
used in this series of papers yields exactly the same finite part as the series of steps carefully developed in [T]. In 
addition, using dimensional regularization to evaluate the singular integrals makes the calculations even more efficient 
as the power-divergent, singular integrals subsequently vanish. We expect that our regularization prescription will 
be most beneficial, if not crucial, for deriving higher-order self force expressions in the gravitational description of 
EMRIs. 

All of the results contained in this paper and in Paper 1 can be derived from one quantity, which is the regular part 
of the master source Su{z'^). The master source can be integrated with the retarded Green's function to generate 
the (scalar) waveform that would be measured by an observer far away, it can be integrated with the regular part of 
the retarded Green's function Dfi{x, z^) to yield the regular part of the field on the worldline and, consequently, the 
regular part of the self force on the SCO. It appears that most, if not all, relevant physical quantities can be derived 
by calculating the master source up to the desired order in e. We also showed how the master source can be calculated 
in perturbation theory using Feynman diagrams and proceeded to demonstrate this with a derivation of the master 
source through fourth order in e. 

We conjecture that a master source S'^(x^) exists for the gravitational description of EMRIs. Preliminary calcula- 
tions suggest that a change of variable from hap to another field Hap can be used to eliminate all cubic self-interaction 
terms in the action that are proportional to □(ft-a^g — (7q^/i^/2), just as in the nonlinear scalar model. It seems likely 
that such a change of variable can be constructed perturbatively to remove higher order interaction terms proportional 
to U{hap — 50/3/1^/2). Therefore, using this field redefinition together with the Detweiler- Whiting decomposition of 
Green's functions into regular and singular parts, we expect that the gravitational perturbations generated by the 
SCO in the new variable Hap may have a representation, at all orders in e, of the form 

Ho.fi{x)^ [ D^pys'{x,x')Sl^'ix') (7.2) 



and expressed in terms of a regular, master source function 5^ (a;^), at least in the Lorenz gauge for trace-reversed 
metric perturbations in a vacuum, background spacetime. If true, one will be able to derive the gravitational waveform 
measured by an observer, the regular part of the metric perturbations H^Jz^^) (and its derivative) on the worldline 
of the SCO, the regular part of the self force, gauge-invariant variables when restricting to conservative dynamics, 
etc., by computing the perturbative expression for essentially one quantity, the master source, to the desired order in 
e. 
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Appendix A: Feynman rules for computing scalar perturbations 

The Feynman rules are adapted (and can be derived) from Paper 1 for the specific purpose of computing scalar 
perturbations. Since we will not be computing an effective action in this paper then we do not need to consider two 
sets of histories for the worldline and field, etc., as outlined in Paper 1. 

The Feynman rules used in this paper are as follows: 
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• For each worldline vertex represented by a dark circle and labeled by c„ write down a factor of {—mcn)- These 
are the worldline vertices. 

• For each curly line write down a factor of D ^{z^, z^ ) = D^ct{z^ , ^^ ) connecting worldline vertices labeled by 
"— " and "+" at proper times r and r', respectively. These are the propagators. 

• Integrate over all proper times. 

• Divide by the appropriate symmetry factor. 

As an example, apply these rules to the diagram on the right in Figure [3] The Feynman rules listed above imply that 
the diagram equals 



2! 



dTdT'dT"D-+{x,z^'){-mc3)D-+{z'',z'' ){-mci)D-+{z^\z^' )(-toci) 



(Al) 



The first factor is the reciprocal of the symmetry factor, which is 2! because there are two ways to connect the ends 
of the two propagators emanating from the C3 vertex to the two ci vertices. The first propagator D '"(x, z'^) in ( Al ) 
is the retarded propagator connecting the C3 vertex at z'^(t) to the field point at x^. The factor of (— mca) is the 
worldline vertex where the three propagators all connect to the worldline. There are two propagators stemming from 
the C3 vertex and connecting from z'^(r) back to the worldline at proper times r' and r" through the two (— ttici) 
worldline vertices. This example shows how one can use the Feynman rules listed in this Appendix to construct any 
such diagram for computing scalar perturbations at any order in perturbation theory. Power counting ( Al ) using the 
scaling laws in Section In] reveals that the diagram on the right in Figure Is] scales as (see (All) 



- (7^)3(7^-2)(m)(7^-2)(m)(7^-2)(TO) 
and is thus a contribution at third order in perturbation theory. 



7^3 



(A2) 



Appendix B: Detweiler- Whiting scheme at all orders in perturbation theory 

In this Appendix, we prove that a regular expression for the master source can be found through any order in 
perturbation theory. We will use dimensional regularization to evaluate the singular integral A, in which case A = 0. 
From the form of the action in (1.7) it is clear that all interactions, by construction from the field redefinition, are 



confined to be on the worldline. Since solutions involving integrals of propagators evaluated at the same point (e.g., 
D{z'^{t), z^^iT))) cannot be generated in a classical theory then every singular contribution to the master source 
contains divergent integrals from the set 



dT2 



■dT,,D,,,{z^\z^-)---D,,t{z^'-\z^-)\ e V)(^''') 



N 



(Bl) 



n=l 



where z^^^ = z^{Ti) for some positive integer i. Let the n'^ member of this set be called J„(2;^i). For example, 5(3) 
in (6.151 contains a J2{z^'^) integral and two Ji(z^i) integrals. 



Consider J„(z''i), any member in the set of (Bl ), and first regularize the t„ integral. In the DW scheme, one writes 
the retarded propagator as 



D,et{x, x') = Ds{x, x') + Dr{x, x') 



(B2) 



so that 



Jn{z^') 



k=l 



dr2---dr„_i( []i^,et(^^^^^'=+^)) I dTn Dn{z^'-\z''-)+Ds{z^'"-\z''-: 

Iniz^"-^) thereby yielding 



From (3.6 1 we know that the t„ integral is just A/(47r) + Ii^{z^^"-^ 

■ n— 3 



Jniz^' 



dT-}-- ■ dr„ 



n ^-t(^ 



MA; ~/^fc + l~ 



fc = l 



dr„ 



Bniz^- 



2^-1)+ 2^^(2^.-2 ^^M„-l) 



(B3) 



Ir{z^-^) (B4) 
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The r„_i integral is then 



(B5) 



The second term involves a proper time integral over Dg and a regular function of t„_i. Since the latter is regular 
we may expand it in a Taylor series for s„_i = t„_i — t„_2 near zero 



/^(z'^-O = luiz^-') + s„_i4(z^'--) + 0(4_i) 



where a dot represents d/dTn-2- Using (3.7) and the above expression, the second term in (B5) becomes 



Since 






dSn-l . " \ (5(Sn-i) = 



for p > 1 then it follows that ( B7 ) equals 



A 



dTn-lDs{z^-\z^'-')lR{z^-^) = —Ir{z^"-') = 

477 



Thus, (B4) is 



J„(Z^1) = JdT2--- dT,,_2 ( n Act(z''Nz'"'+0) j dT^-1 Dr{ 



z^-\z'^-')Ir{z^-') 



The remaining proper time integrals are evaluated in like manner. In particular, by induction it follows that 



j„(z^^; 



j dT2---dTj\{ Dr{z''^ 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 



and the regular master source is constructed from integrals in the set appearing in the right side of (Bl) with D^ct 
replaced by the regular part. Dr. Hence, one can always find the regular part of the master source Sr(z^^^) to any 
order in perturbation theory. 

Now for the self force. Since we have shown that the regular part of the master source can always be constructed 
at any order in perturbation theory then the field is simply 



i,{x) = 4'r{x) = / dr'DRix.z^^ )SR{z^'') 



(B12) 

which is regular on the worldline when a;^ = ^^{.t)- The self force also depends on the derivative of ^ on the worldline. 
Using calculations presented in Paper 1 it is straightforward to show that 



47r z 



Rewriting the worldline equations of motion in (1.10) as 



ma^' = -m{a^' + P^''V^)C(V'(z^)) = F''(t) 



it follows from (B12) and (B14) that the self force on the small compact object is 

F^{t) = ~m{a^ + P^'^VjC(^fl(z'^)) = F]^(t) 



(B13) 
(B14) 

(B15) 

(B16) 



Therefore, the Detweiler- Whiting scheme for calculating the self force is valid at any order in perturbation theory, at 
least in this class of nonlinear scalar theories. Of course, V'fi(z^) must be computed perturbatively and we performed 
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the calculation explicitly in this paper through third order in e and the master source through fourth order. Indeed, 
the EFT approach provides a systematic way of computing ipn to any order in perturbation theory using Feynman 
diagrams. 
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